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Saturation  threshold  of  a  probe  pulse  in  an  ultrafast  electronic-resonance-enhanced  CARS  (ERE- 
CARS)  configuration  is  calculated.  We  demonstrate  that  while  underdamping  condition  is  a  suffi¬ 
cient  condition  for  saturation  of  ERE-CARS  with  the  long-pulse  excitations,  a  transient-gain  must 
be  achieved  to  saturate  ERE-CARS  signal  for  ultrafast  probe  regime.  We  have  identified  that  the 
area  under  the  probe  pulse  can  be  used  as  a  definitive  parameter  to  determine  the  criterion  for 
saturation  threshold  for  ultrafast  ERE-CARS.  From  a  simplified  analytical  solution  and  a  detailed 
numerical  calculation  based  on  density-matrix  equations,  the  saturation  threshold  of  ERE-CARS  is 
compared  for  a  wide  range  of  probe-pulse  durations  from  10-ns  to  10-fs  regime.  The  theory  explains 
both  qualitatively  and  quantitatively  the  saturation  thresholds  of  resonant  transitions,  and  also 
gives  a  predictive  capability  for  other  pulse  duration  regimes.  The  presented  criterion  for  saturation 
threshold  will  be  useful  in  establishing  the  design  parameters  for  ultrafast  ERE-CARS. 

PACS  numbers:  42.65.Dr,  42.65.-k,  82.53.Kp 


I.  INTRODUCTION 

Femtosecond  (fs)-laser-based  coherent  anti-Stokes  Ra¬ 
man  spectroscopy  (CARS)  has  been  used  extensively  in 
time-resolved  nonlinear  spectroscopy  research  in  recent 
years  [1].  A  variant  of  CARS,  a  doubly- resonant  four- 
wave  mixing  scheme,  where  a  probe  pulse  resonant  to 
the  electronic  transition  scatters  off  the  Raman  coher¬ 
ence  generated  by  the  Raman-resonant  pump  and  Stokes 
pulses,  is  known  as  electronic-resonance-enhanced  CARS 
(ERE-CARS)  [2] .  This  technique  is  capable  of  obtaining 
a  few  orders  of  magnitude  enhancement  of  CARS  signal 
as  compared  to  traditional  CARS  with  a  non-resonant 
probe  [3-5].  Also,  since  ERE-CARS  technique  is  inher¬ 
ently  quenching  independent  [6,  7],  this  technique  has 
been  employed  very  successfully  using  nanosecond  (ns) 
lasers  for  measuring  concentration  of  minor  species  such 
as  NO  [2,  6,  8-12],  which  is  a  tracer  pollutant  in  combus¬ 
tion  processes  [13].  However,  since  the  resonant  probe 
in  ERE-CARS  configuration  interacts  strongly  with  the 
molecules,  the  saturation  limit  of  the  probe  intensity 
is  lower  than  that  in  conventional  CARS  [14].  It  has 
been  observed  that  to  maximize  the  signal-to-noise  ratio 
(SNR),  the  ns  ERE-CARS  setup  were  operating  near  sat¬ 
uration  limit  of  the  probe  intensity  [16].  In  such  spectro¬ 
scopic  techniques,  while  it  is  important  to  obtain  spec¬ 
trum  of  the  intended  molecule  in  the  generated  ERE- 
CARS  signal  with  acceptable  SNR,  saturation  can  al¬ 
ter  the  spectral  signature  (such  as  Stark  shift  of  tran¬ 
sitions  [15]),  requiring  additional  model  corrections  like 
those  employed  in  saturation  of  fluorescence  [17],  four- 
wave  mixing  [18,  19]  and  standard  off-resonant  CARS 
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[20].  With  the  ultrafast-laser  excitations,  since  the  oper¬ 
ating  intensities  of  the  lasers  are  generally  much  higher 
than  those  used  in  ns  ERE-CARS  experiments  and  the 
duration  of  exciting  pulses  are  much  shorter  than  the 
molecular  decay  and  dephasing  timescales,  it  is  very  im¬ 
portant  to  understand  the  ultrafast  saturation  dynamics 
of  the  ERE-CARS. 

Recently,  we  demonstrated  that  the  physics  of  ultra¬ 
fast  saturation  of  elementary  optical  processes  (such  as 
absorption  and  fluorescence)  differs  significantly  for  dif¬ 
ferent  durations  of  the  pumping  laser  pulses  [21,  22]. 
It  was  shown  that  in  an  underdamping  regime,  i.e. ,  if 
the  Rabi  frequency  associated  with  a  long-pulse  or  CW 
laser  is  greater  than  the  effective  decay  rate  of  the  sys¬ 
tem,  absorption  and  fluorescence  saturation  can  occur 
[23,  24].  However,  for  short  pulses  when  the  pulse  dura¬ 
tion  of  the  pump  is  shorter  than  the  effective  decoherence 
timescales  of  molecule,  the  saturation  of  ultrafast  absorp¬ 
tion  occurs  only  when  the  pulse  drives  the  molecule  to  a 
transient-gain  regime.  Such  differences  in  the  saturation 
of  polarization  spectroscopy  with  long  and  short  pulses 
have  been  reported  earlier  [25,  26].  The  ultrafast  satu¬ 
ration  limit  could  be  predicted  quantitatively  employing 
the  pulse  area  [27]  as  the  deterministic  parameter,  and 
the  same  argument  is  extended  to  saturation  of  the  Ra¬ 
man  transition  [22].  Note  that  the  saturation  threshold 
of  probe-pulse  intensity  in  a  traditional  CARS  is  much 
higher  than  the  ERE-CARS  setup  because  of  far  weaker 
off-resonant  interaction  of  the  probe  with  the  molecular 
transition  [21].  Although  saturation  of  ERE-CARS  with 
long-pulse  excitation  has  been  extensively  studied  [7,  14] , 
the  saturation  of  the  ultrafast  ERE-CARS  is  never  dis¬ 
cussed. 

In  this  paper,  we  generalize  study  of  ultrafast  satu¬ 
ration  of  emission  for  quantitative  determination  of  the 
threshold  of  saturation  of  ultrafast  ERE-CARS.  Satura¬ 
tion  of  the  ultrafast-CARS  signal  with  an  off-resonant 
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probe  setup  was  discussed  in  Ref.  [21].  Here,  we  obtain 
a  quantitative  saturation  criterion  for  the  probe  inten¬ 
sity  to  determine  the  threshold  of  saturation  of  ultra¬ 
fast  ERE-CARS  using  the  area  under  the  probe  pulse. 
We  compare  the  calculated  threshold  of  saturation  for 
a  range  of  short-  and  long-pulse-duration  probe  pulses 
to  compare  and  fundamentally  distinguish  the  satura¬ 
tion  processes  for  those  durations.  We  show  that  the 
pulse-area-based  criterion  for  saturation  of  ERE-CARS 
is  valid  for  shorter  pulses  and  is  accurate  for  the  delayed 
probe  conditions.  Our  study  is  valid  for  ultrashort  pulses 
where  rotating-wave  approximation  (RWA)  is  still  valid, 
i.e. ,  the  Rabi  frequency  is  much  smaller  than  the  optical 
frequency  [28].  Also,  the  propagation  effect  is  not  con¬ 
sidered.  Note  that  the  ultrafast  saturation  reported  here 
is  on  a  time  scale  of  the  vibrational  period  of  molecules 
[29]  and  is  orders  of  magnitude  faster  than  state-of-art 
electronic  switches.  This  study  provides  predictive  ca¬ 
pability  for  the  intensity  threshold  to  avoid  ERE-CARS 
saturation,  which  is  very  important  for  establishing  de¬ 
sign  parameters  needed  for  ultrafast  spectroscopic  mea¬ 
surements.  Other  interesting  applications  may  include 
saturation-based  fs  switches  and  controllable  wave-packet 
dynamics  [30],  above-saturation-threshold  spectroscopy 
[31] ,  and  enhanced  spatial  and  spectral  resolution  via  sat¬ 
urated  CARS  [32]. 

The  paper  is  organized  as  follows:  in  Sec  II,  we  dis¬ 
cuss  the  model  and  derive  a  density  matrix  equation  to 
calculate  ERE-CARS  polarization.  In  Sec.  Ill,  we  ob¬ 
tain  approximate  analytical  and  numerical  solutions  to 
understand  the  ultrafast  dynamics  of  the  molecule  under 
the  influence  of  the  strong  probe  laser.  The  ultrafast  sat¬ 
uration  of  signal  in  a  hybrid  ERE-CARS  configuration  is 
presented  and  effect  of  probe  delay  on  saturation  is  dis¬ 
cussed  in  Sec.  IV.  Comparison  of  saturation  criterion  is 
obtained  numerically  and  compared  with  analytical  re¬ 
sults  for  a  wide  range  of  probe  pulse  durations  in  Sec.  V. 
The  results  are  summarized  in  Sec.  VI. 


II.  THE  MODEL  CALCULATION  FOR 
ERE-CARS  POLARIZATION 

We  consider  a  six-level  system  for  studying  the  satu¬ 
ration  of  ERE-CARS  from  ns  to  fs  excitation  regimes,  as 
depicted  in  Fig.  1.  The  system  parameters  used  for  the 
calculation  are  kept  close  to  that  of  an  NO  molecule.  The 
energy  levels  |a)  and  la')  could  correspond  to  a  rovibra- 
tional  state  in  the  excited  electronic  state  A2E+(u  =  0). 
The  energy  levels  |6)  and  |6')  (|c)  and  |c'))  are  the  rovi- 
brational  states  in  the  v  =  1  (v  =  0)  vibrational  mani¬ 
folds  in  the  ground  electronic  state  X2n.  The  primed 
states  are  representative  states  of  the  rotational  manifold 
within  the  same  vibrational  states  as  the  corresponding 
unprimed  states,  and  they  are  used  as  bath  levels  for 
accurately  taking  into  account  the  various  decays  such 
as  rotational-energy  transfer  (RET),  vibrational-energy 
transfer  (VET),  electronic  decay,  and  quenching  rates  of 


FIG.  1.  (Color  online)  Six-level  model  used  for  studying 
ERE-CARS  saturation  of  NO.  The  rotational  manifolds  in 
the  ground  (excited)  vibrational  state  in  the  ground  electronic 
state  X2n  are  represented  by  |c)  and  \c'),  and  rovibrational 
states  in  the  excited  electronic  state  A2E+  are  represented  by 
|a)  and  |a').  The  pump,  Stokes,  probe,  and  generated  CARS 
pulses  that  are  represented  by  E i,  E2,  E3  and  E4,  respec¬ 
tively,  couple  the  unprimed  states  |a),  |6),  and  |c).  The  prime 
states  are  considered  to  be  the  bath  levels  to  account  for  all  of 
the  decays  (IV7)  in  the  system  and  the  coherence  dephasing 
is  represented  by  7 ij. 


the  NO  molecule.  The  population-decay  and  coherence¬ 
dephasing  timescales  vary  in  the  ns  to  ps  range;  hence, 
the  corresponding  rates  are  important  in  attaining  ac¬ 
curate  ERE-CARS  results  in  the  ns-  and  ps-excitation 
regimes.  However,  the  decays  and  dephasing  rates  play 
very  insignificant  role  in  time-resolved  signal  detection  in 
the  fs-excitation  regime.  For  a  detailed  discussion  on  the 
decays  and  dephasing  parameters  used  in  the  modeling 
of  the  NO  molecule,  see  Ref.  [7].  Furthermore,  note  that 
a  short  pulse  may  couple  multiple  rotational  states  in 
each  vibrational  manifold  simultaneously  because  of  the 
broad  bandwidth  associated  with  the  pulse  that  affects 
the  saturation  dynamics  of  the  CARS  signal  [21].  How¬ 
ever,  since  the  goal  of  this  work  is  to  obtain  an  analytical 
understanding  of  the  saturation  threshold  of  ERE-CARS 
for  pulse  durations  ranging  from  ns  to  fs,  we  limit  our 
model  with  pump,  Stokes,  and  probe  coupling  to  only 
three  of  the  molecular  levels. 

Three  input  pulses  Ei  (i  —*■  1,2,3)  at  frequencies  f, 
generate  the  CARS  signal  E4  at  a  frequency  1/4  =  v\  — 
V2+V3-  The  pump  (E\)  and  Stokes  ( Ei )  pulses  are  in  two- 
photon  resonance  with  the  Raman  transition  |  b)  o  |c)  to 
generate  the  Raman  coherence  pbc ■  The  probe  E3  couples 
to  the  electronic  transition  via  |a)  -o-  |6)  resonantly  or 
near  resonantly  with  a  detuning  of  A.  All  of  the  input 
pulses  are  given  by 

Ei(t)  =  eiEi(t)e~Wit  +  c.c.,  (1) 

with  Gauusian  envelope  functions 

Ei(t)  =  El0e~^2^ . 


(2) 


3 


Here,  E,q ,  e*,  f*,  and  rt  are  the  peak  electric  field,  polar¬ 
ization,  peak  position  of  the  pulse,  and  duration  of  the 
pulse  of  the  applied  fields,  respectively  (z  — ►  1, 2, 3).  All 
of  the  fields  are  assumed  to  be  linearly  polarized  and  are 
parallel  to  each  other  and  parallel  to  the  induced-dipole 
moments  corresponding  to  their  respective  transitions. 
The  total  Hamiltonian  is  given  by 

H  =  h(u;  ac-A-aa  A  LObc-Abb  A  (0a'cAa'a/  A  bJb'cAb'b1)  (3) 
-  H  \p.\e~ivit Aac  A  (H2e-iiy2t  A  n3e~iV3t)Aab  A  H.c]  , 

Here,  huiac  is  the  energy  separation  between  the  states  a 
and  the  ground  state  c  (a  — >  a,  b,  a',  b').  The  projection 
operators  Aap  =  |a)(/3|  correspond  to  population  (coher¬ 
ence)  for  a  =  ^  (a  /  /3),  where  a,/3  — >  a,b,c,a',b',c'. 
Note  that  the  prime  states  are  assumed  to  be  degenerate 
with  the  corresponding  unprimed  states  for  simplicity. 
The  Rabi  frequencies  corresponding  to  the  incident  fields 
are  denoted  as 

Hi  —  dac •  e i  Ei(t)/H 
=  dab-e2  E2(t)/h 

H3  =  dab-e 3  Ez(t)/h.  (4) 

Here  dab  ( dac  )  is  the  transition  matrix  element  that  in¬ 
cludes  the  induced  dipole  moment,  Frank-Condon  factors 
and  the  rotational-overlapping  function  between  the  rovi- 
brational  levels  |a)  —>  |6)  (|a)  —>  |c)). 

The  Liouville  equation  to  describe  the  dynamics  of  the 
molecule-laser  interaction  for  the  model  system  can  be 
obtained  as 

%  =  -^[H,p(t)]+Cp(t),  (5) 

where  C  is  the  Liouville  operator  consisting  of  the  the 
phenomenological  decay  (dephasing)  rates  of  the  molec¬ 
ular  state  population  (coherences).  Substituting  Eq.  (4) 
into  Eq.  (5)  and  writing  explicitly  all  of  the  phenomeno¬ 
logical  decay  and  depliasing  rates  of  the  system,  the  equa¬ 
tions  for  the  density-matrix  elements  can  be  written  as 

Paa  —  Paa  A  Ta/a  Pa'a'  A  [z  Hie  Pea 
Az  (H2e  ll'2t  A  H3e  1V3t)pba  A  c.c.] , 

Pab  ( 7a6  A  ZC dab')  Pab  A  Z  Hie  pbc 

Az  (H2e  lV2t  A  Hse  lV3t)ipbb  —  Paa), 

Pac  —  {'lac  A  IbJac)  Pac  A  Z  (H2e  2  A  H3e  3  )  pbc 
AZ  Hie  1  ( pec  Paa): 

Pbb  =  — r bb>  Pbb  ATm  pb’b' 

-  [z  (H2e_“'2*  A  H3e_W3t)pa6  A  c.c.] ,  (6) 

Pbc  =  -(76c  -  Z u>bc)  Pbc  A  z  (H2eW2t  A  ^leW3t)pac 

-tHie-^Vfco, 

Pcc  r Cc '  Pcc  A  rc'c  Pc'c'  [z  Hie  Pea  A  C.C.]  , 

Pa'a'  —  ha'  Pa'a'  A  Taa'  paa: 

Pb'b'  —  h b'b  Pb'b'  A  Ta5/  Paa  A  T a'b'  Pa'a'  A  Tbfr/  Pbb: 

Pc'c'  r e'e  Pc'c'  A  F ac'  Paa  A  T a'c'  Pa'a'  A  F Cc'  Pcc  • 


Note  that  the  density-matrix  elements  pap  represent  the 
population  (coherence)  for  a  =  /3  (a^/3),  and  =  ppa. 
The  population  decay  rates  associated  with  the  excited 
electronic  states  Ta  and  Fa/  are  the  sum  of  the  contri¬ 
butions  from  spontaneous  decay,  RET  and  collisional 
quenching  rates.  The  decay  rate  Faa/(=  10  x  Fa'a) 
is  considered  to  be  equal  to  the  RET  rate.  The  de¬ 
cay  rate  corresponding  to  the  electronic  transitions  is 
7 aj  =  (ro  +Tjj>)/2  +  Jaj  U  b,c)  and  that  correspond¬ 
ing  to  the  Raman  transition  7 bc  is  primarily  the  RET 
rate.  For  details  of  the  decay  model  considered  here,  see 
Ref.  [7], 

To  separate  the  molecular  polarizations  oscillating  at 
different  frequencies  and  to  remove  the  fast  oscillating 
terms  in  the  density-matrix  equations,  we  used  the  fol¬ 
lowing  canonical  transformations: 

Pac{t)  =  (Tacit) e~Wlt  A  Pac{t)e~lVit , 

Pab{t)  =  aab{t)e~tV2t  A  r]abit)e~ll'3t ,  (7) 

Pbc(t)  =  Pbc(t)e-«Vl-^\ 

Here  zq  is  the  frequency  of  the  generated  ERE-CARS 
signal.  To  obtain  the  equations  in  transformed  frame, 
we  dropped  the  fast  terms  oscillating  at  e±i<5t  with  the 
assumption  that  <5  ^  f  .  Note  that  16  coupled  differ¬ 
ential  equations  are  obtained  after  substituting  Eqs.  (7) 
into  Eqs.  (6)  that  describe  the  dynamical  evolution  of 
the  the  molecular  states  and  coherences. 


III.  ULTRAFAST  EVOLUTION  IN  ERE-CARS 
CONFIGURATION 

To  obtain  an  analytical  understanding  of  the  evolution 
of  ERE-CARS  polarization,  let  us  evaluate  the  coherence 
Tjac  that  oscillates  at  zq,  representing  the  instantaneous 
CARS  polarization.  The  equation  of  dynamics  for  the 
matrix  element  r/ac  after  applying  the  transformations 

(7)  to  Eqs.  (6)  is 

d>la^  =  _  (7oc  +  iA)  pa c(t)  A  zH3(t)pfec(t).  (8) 
ot 

Here,  7 ac  =  \ila  +  lc)  +^aF  is  the  total  coherence  decay 
associated  with  the  |a)  |c)  transition.  Expanding  Eq. 

(8)  in  a  Taylor  series  for  H3  — >  0,  a  formal  solution  of  the 
coherence  to  the  first  order  in  H3  can  be  obtained  using 
an  integrating  factor  [33] 


a(t)  =  Exp 


(7 ac  A  iA)dt' 


UO 


(9) 


which  leads  to  the  formal  solution  for  coherence  pac  to 
the  first  order  in  the  probe  [34] 

V&Ht)  =  ze-^-+iA)‘  f  H3(^')/5i0c)(^,)e(7“•=+^A)t'd^,. 

J 1 30—  T3 

(10) 
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If  the  probe  delay  is  much  longer  than  the  pulse  durations 
of  the  pump  and  Stokes,  we  assume  that  the  probe  delay 
is  much  longer  than  t\  and  £2!  then  we  can  assume 
that  the  zeroth-order  Raman  coherence  (in  absence  of  the 
probe)  remains  constant  within  the  duration  of  the  probe, 
i.e., 


Pbctt)  =  Pbc  (say)  (11) 


is  a  constant  within  the  limit  £3  —  13  <  t  <C  q//1 .  For  the 
Gaussian  probe  appearing  within  this  time  window,  the 
coherence  in  the  ERE-CARS  transition  can  be  obtained 
(with  A  =  0)  as 


??W(i)  =  (i/2)V^e^+^Ts/2)2  Pbcn^ne-^ 

't  ^30  7ocT3 


Erf  (l  +  +  Erf 


73 


(12) 


Note  that  the  above  coherence  is  purely  imaginary  be¬ 
cause  we  have  considered  A  =  0.  At  a  longer  time 
limit,  t  — >  7//1,  the  error  functions  remain  unchanged 
but  the  coherence  rjac(t)  decays  to  become  zero  at  a  time 
t  7 ac-  However,  for  an  intense  probe,  the  above  result 
will  become  invalid  because  the  intense  probe  can  perturb 
the  ground-state  coherence  pbc,  as  will  be  discussed  later 
with  reference  to  Fig.  2(c)  obtained  from  the  numerical 
solutions. 

Under  the  same  condition  as  above  and  assuming  the 
limiting  case  of  7 ac  0  and  (£3  —  73)  — >  —00,  the  integral 
in  Eq.  (8)  reduces  to  the  transient  pulse  area  of  the  probe 
pulse  63 (t);  thus, 


VaJ  =  (  7  )  Pbc03(t), 
where  63  is  defined  as  [22] 


(13) 


However,  note  that  the  assumption  (11)  that  is  used  to 
obtain  the  above  solution  clearly  indicates  that  the  pulse- 
area-based  solution  may  not  be  valid  for  probe  durations 
t3  >  7 be  ■  In  the  later  part  of  the  paper,  we  show  from 
our  numerical  calculations  that  above  saturation  condi¬ 
tion  is  indeed  not  valid  if  the  pulses  are  longer  than  the 
effective  decay  time  scale.  It  has  been  shown  earlier  that 
with  ultrafast-laser  excitations,  the  absorption  and  flu¬ 
orescence  in  a  two-level  system  saturate  when  the  area 
under  the  driving  pulse  reaches  7r/2  and  7 r,  respectively 
[22]. 
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FIG.  2.  (Color  online)  (a)  Population  in  excited  electronic 
state  |a),  (b)  excited  vibrational  state  |fc),  (c)  ground-state 
(Raman)  coherence,  and  (d)  instantaneous  polarization  in  the 
ERE-CARS  transition  for  different  probe  intensities.  The  pa¬ 
rameters  used  for  the  plot  are:  the  peak  of  the  pump  and 
Stokes  appears  at  400  fs,  and  the  probe  delay  is  maintained 
at  200  fs. 


03(t)=  f  2 n3(t')dt'.  (14) 

J  —00 

From  Eq.  (13)  it  is  clear  that  the  instantaneous  ERE- 
CARS  polarization  is  proportional  to  the  pulse  area  of 
the  probe  field  under  the  weak-probe  condition.  For  the 
Gaussian  probe  pulse  under  consideration,  the  transient 
polarization  r^c  can  be  rewritten  as 

VaJ  =  iV^Pbc^3  ORJ-  (15) 

Here,  H30  is  the  peak  Rabi  frequency  over  the  duration 
of  the  probe  pulse.  Thus,  the  coherence  r/aJ  is  propor¬ 
tional  to  73.  Even  though  the  above  result  is  obtained 
in  weak-probe  condition,  the  result  depicts  the  under¬ 
lying  physics  for  existence  of  a  threshold  criterion  for 
saturation  of  ERE-CARS  signal.  Since  the  value  of  any 
coherence  can  reach  a  maximum  value  of  only  0.5,  the 
threshold  Rabi  frequency  of  the  probe  for  saturation  of 
is  inversely  proportional  to  the  duration  of  the  pulse. 


In  real  ERE-CARS  experiment,  the  probe  must  be 
much  stronger  if  an  acceptable  SNR  is  to  be  achieved. 
Hence,  to  obtain  a  solution  for  ERE-CARS  polarization 
under  strong  probe  conditions,  we  obtain  a  numerical 
solution  to  study  the  ERE-CARS  dynamics.  We  con¬ 
sider  the  dynamics  of  the  molecular-state  population  ( paa 
and  pbb)  and  Raman  coherence  (pbc)  for  ultrafast  pump, 
Stokes,  and  probe  pulses  and  also  determine  how  the 
instantaneous  ERE-CARS  polarization  rjac  (and,  hence, 
the  ERE-CARS  signal)  evolves  with  the  ultrafast  pulses. 
The  coupled  density-matrix  equations  are  solved  numer¬ 
ically  using  a  fifth-order  Runge-Kutta  method,  and  the 
results  are  presented  in  Fig.  2.  While  keeping  the  values 
of  the  peak  intensities  of  the  pump  and  Stokes  pulses  un¬ 
changed  at  J10  =  ho  =  1010  W/m2,  we  varied  the  peak 
intensity  of  the  probe  field  I30  parametrically,  ranging 
between  1011  and  1018  W/m2.  For  all  of  the  numerical 
results,  we  considered  the  pump,  Stokes,  and  probe  du¬ 
rations  to  be  Ti  =  100  fs.  The  peak  of  the  pump  and 
Stokes  appears  at  t\  =  t2  =  400  fs.  As  in  a  typical 
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fs-CARS  experiment,  we  consider  the  probe  delay  to  be 
200  fs  and,  hence,  t3  =  600  fs.  Note  that  assuming  the 
beam  diameter  to  be  100  p m,  the  peak  intensity  of  the 
probe  I3 o  =  1011  W/m2  correspond  to  the  laser  energies 
of  1.4  nanoJoule  (nJ),  140  picoJoule  (pj)  and  14  pj  at 
probe  durations  r3  =  1  ps,  100  fs  and  10  fs,  respectively. 
The  maximum  value  of  /3 o  =  1018  W/m2  corresponds 
to  laser  energies  of  14  miliJoule  (mJ),  1.4  mJ  and  140 
/iJ  for  r3  =  1  ps,  100  fs  and  10  fs,  respectively.  Such 
laser  energies  are  readily  available  in  commercial  high- 
power  lasers.  For  example,  a  Ti:sapphire  laser  operating 
at  800  nm  at  6  mJ  per  pulse  laser  energy  generating  100 
fs  pulses  can  produce  about  20  pj  at  236  nm  using  a 
tunable  optical  parametric  amplifier. 

Figure  2(a)  shows  that  for  the  weaker  probe  /3o  <  1014 
W/m2,  the  population  in  the  excited  electronic  state  paa 
reaches  a  maximum  at  ~600  fs,  i.e. ,  when  the  peak  of 
the  probe  pulse  appears.  However,  since  the  duration 
of  the  pulse  is  very  short  compared  to  the  molecular  re¬ 
sponse  time,  the  excited-state  population  does  not  follow 
the  probe  pulse  during  the  tail  end  of  the  pulse;  later,  paa 
decays  with  a  characteristic  decay  rate  Ta  of  the  state  |a). 
For  higher  probe  intensities,  the  transition  saturates  even 
before  the  peak  of  the  probe  pulse  appears  at  f3;  hence, 
the  maximum  achievable  value  of  paa  is  reduced.  For 
ho  >  1014  W/m2,  the  Rabi  oscillation  is  observed  in  paa 
after  the  onset  of  saturation.  Once  the  probe  pulse  leaves 
the  molecule,  the  population  decays  at  the  rate  Ta.  Simi¬ 
larly,  the  population  pbb  in  the  excited  rovibrational  state 
|6)  (in  the  ground  electronic  state)  in  Fig.  2(b)  exhibits 
a  monotonic  increase  initially  when  the  pump  and  Stokes 
pulses  are  present.  It  reaches  a  quasi-steady-state  for  the 
time  interval  t3  <  t  <C  7//1  in  the  absence  of  the  probe 
(not  shown  here).  However,  once  the  probe  is  turned  on 
it  resonantly  couples  the  transition  |a)  -A  |c)  by  orders- 
of-magnitude  stronger  coupling  than  the  Raman  coupling 
by  the  pump  and  Stokes  pulses  with  \b)  A  |c).  Hence, 
the  probe  takes  away  molecular  population  from  state  |  b) 
to  populate  state  |a).  Thus,  pbb  remains  low  compared  to 
Paai  even  for  weaker  probe  intensities.  However,  once  the 
probe  transition  is  saturated  and  after  the  onset  of  Rabi 
oscillation,  the  population  oscillates  between  |a)  and  |6), 
but  their  values  remain  out  of  phase  by  7 r.  Note  that 
for  the  delayed  resonant  probe  considered  here,  the  pop¬ 
ulation  in  the  probe  transition  oscillates  as  it  does  in  a 
resonant  two- level  system  [22]. 

The  ground-state  coherence  \pbc\  is  also  plotted  in 
Fig.  2(c),  which  shows  behavior  similar  to  that  of  pbb', 
however,  they  are  out  of  phase  by  7 r.  For  stronger 
probe  pulses,  pbc  shows  Raman-saturation-like  oscilla¬ 
tions,  which  were  described  by  the  two-photon  pulse  area 
approach  [21,  22].  However,  note  that  the  origin  of  the 
ultrafast  oscillations  in  those  references  was  because  of 
saturation  of  the  Raman  transition  by  strong  pump  and 
Stokes  pulses,  unlike  in  the  current  study  where  the  am¬ 
plitude  of  Raman  coherence  is  well  below  its  maximum 
value  of  0.5.  Here,  the  resonant  probe  pulse  causes  os¬ 
cillation  of  the  population  in  the  state  |6),  resulting  in 
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FIG.  3.  (Color  online)  Saturation  of  ERE-CARS  signal  inten¬ 
sity  for  different  pulse  durations  in  hybrid  ERE-CARS  con¬ 
figuration.  Pump-  and  Stokes-pulse  durations  are  maintained 
at  n  =  T2  =  100  fs,  and  only  the  probe-pulse  duration  r3 
is  varied,  (a)  ERE-CARS  signal  Sere  is  plotted  as  a  func¬ 
tion  of  peak  intensity  of  the  probe  /30,  and  (b)  pulse  area  ©3 
corresponding  to  ho- 


oscillation  in  pbc-  In  Fig.  2(d),  we  present  the  coherence 
\r]ac(t)\  that  depicts  the  generated  instantaneous  ERE- 
CARS  polarization  or  the  square  root  of  the  transient 
ERE-CARS  signal.  For  a  weak  probe,  the  coherence 
77 ac  shows  an  asymmetric  structure  in  time,  reaching  the 
maximum  when  the  peak  of  the  probe  pulse  appears.  In 
the  next  section,  we  consider  saturation  of  the  integrated 
ERE-CARS  signal. 


IV.  ULTRAFAST  SATURATION  OF  ERE-CARS 

To  investigate  saturation  of  ERE-CARS  with  a  strong 
probe,  if  we  assume  that  the  excitation  of  paa  in  the 
probe  transition  |  b)  A  |a)  equivalent  to  an  isolated  two- 
level  transition  [22]  and  that  the  ERE-CARS  signal  is 
proportionally  dependent  on  paa,  we  may  predict  that 
the  ERE-CARS  signal  will  saturate  once  total  pulse  area 
reaches  ir,  i.e., 

03  =  lim  d3(t)  =  2n30T3V/7r  >  7 r.  (16) 

t—>  OO 

This  relation  shows  that  the  saturation  threshold  of  Rabi 
frequency  is  inversely  proportional  to  the  pulse  duration; 
and,  hence,  the  threshold  intensity  of  the  probe  for  sat¬ 
uration  of  ERE-CARS  signal  is  quadratically  dependent 
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on  the  inverse  of  the  pulse  duration,  i.e. , 


jth  _ 

i30  — 


23|dgc|2  -2 
h2Trceo  3 


(17) 


where,  c  and  eo  are  the  speed  of  light  and  the  permitivity 
in  vacuum,  respectively.  Note  that  since  the  ERE-CARS 
configuration  is  more  complex  and  the  signal  is  depen¬ 
dent  on  the  Raman  coherence,  we  numerically  investigate 
the  saturation  of  the  ERE-CARS  signal  for  ultrafast  ex¬ 
citation. 

In  Fig.  3,  we  present  the  numerical  results  to  verify  the 
above  prediction.  We  consider  three  different  probe-pulse 
durations  (73  =  1  ps,  100  fs  and  10  fs)  while  maintaining 
the  durations  of  the  pump  and  Stokes  pulses  at  t\  =  72  = 
100  fs.  Typically,  the  CARS  configurations  with  probe 
durations  different  from  those  of  the  pump  and  Stokes 
pulses  are  called  a  hybrid-CARS  configurations  [35-37], 
and  the  current  configuration  may  be  called  as  hybrid- 
ERE-CARS  configuration.  All  parameters  are  assumed 
to  have  the  same  value  as  that  in  Fig.  2.  In  Fig.  3(a), 
we  present  the  numerical  solution  of  the  integral  of  the 
square  of  the  coherence  \r]ac\2  that  represents  integrated 
ERE-CARS  signal 


Sere  oc 


\r]ac{t)\2  dt 


(18) 


as  a  function  of  the  probe  intensity  /30.  Clearly,  the 
threshold  intensity  of  the  probe  for  saturation  of  the 
ERE-CARS  signal  increases  by  two  orders  of  magnitude 
when  the  pulse  duration  of  the  probe  decreases  by  one 
order  of  magnitude.  This  dependence  is  precisely  cap¬ 
tured  in  Eq.  (15),  even  though  Eq.  (15)  was  obtained  in 
the  weak-probe-field  limit. 

Furthermore,  Fig.  3(b)  is  a  numerical  plot  of  the  pulse 
area  corresponding  to  the  probe  intensities  used  in  the 
x-axis  of  Fig.  3(a).  Comparing  the  Figs.  3(a)  and  3(b), 
the  ERE-CARS  saturation  threshold  is  reached  once  the 
pulse  area  corresponding  to  the  probe  area  is  larger  than 
7T.  Note  that  Figs.  3(a)  and  3(b)  are  plotted  on  a  log- 
log  scale.  The  saturation  criterion  agrees  well  with  the 
10-fs  and  100-fs  probe  pulses.  However,  the  03-based 
criterion  slightly  underpredicts  the  threshold  limit  of  the 
Raman  saturation  for  the  ps  probe  by  a  small  fraction 
of  the  actual  Raman  threshold.  The  reason  for  this  is 
that  the  front  end  of  the  ps  probe  significantly  overlaps 
with  the  pump  and  Stokes  pulses,  leading  to  population 
recycling  between  states  \b)  and  |c)  during  the  pulse  ex¬ 
citation.  Hence,  the  03-based  criterion  (17)  is  not  exact 
for  longer  pulses,  which  was  also  mentioned  in  the  discus¬ 
sion  following  Eq.  (15).  However,  it  lends  a  predictive 
capability  to  determine  the  saturation  threshold  of  ultra¬ 
fast  ERE-CARS  within  an  order  of  magnitude  not  only 
for  the  fs  probe  but  also  for  a  range  of  pulse  durations. 

Note  that  the  pulse-area-based  saturation  criterion  is 
derived  under  the  condition  that  the  probe  pulse  is  de¬ 
layed  with  respect  to  the  pump  and  Stokes.  To  under¬ 
stand  how  the  saturation  criterion  obtained  from  ©3  may 
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FIG.  4.  (Color  online)  Saturation  of  ERE-CARS  signal  in¬ 
tensity  for  different  probe  delays. 


be  affected  by  the  probe  delay,  we  plotted  the  ERE- 
CARS  signal  in  Fig.  4  as  a  function  of  I30  for  differ¬ 
ent  probe  delays.  We  consider  the  duration  of  the  probe 
pulse  to  be  73  =  100  fs,  and  all  other  parameters  are 
maintained  at  the  same  value  as  in  Fig.  3.  The  plots 
in  Fig.  4  show  that  if  the  probe-pulse  delay  is  zero  or 
small  {t delay  <  T3),  the  saturation  of  the  ERE-CARS  sig¬ 
nal  occurs  at  a  much  larger  value  of  /30  than  that  for 
rd  ~  200  or  300  fs.  For  100-fs  pulse  duration,  03  =  7r 
for  J30  «  5  x  1014  W/m2.  Thus,  the  03-based  criterion 
underpredicts  the  saturation  threshold  for  small  or  no 
probe  delays.  This  discrepancy  may  be  caused  by  sig¬ 
nificant  repopulation  of  state  |6)  in  the  presence  of  the 
pump  and  Stokes  pulses.  For  a  probe  delay  of  ~  200  fs, 
the  saturation-threshold  intensity  follows  the  pulse-area- 
based  criterion  of  ©3  >  7 r.  Once  again,  the  ©3-based 
criterion  underpredicts  the  threshold  intensity  for  probe 
delays  in  excess  of  400  fs.  The  saturation  threshold  con¬ 
tinue  to  increase  for  even  larger  probe  delays.  Such  a 
discrepancy  could  be  caused  by  a  decrease  in  the  pop¬ 
ulation  in  | b)  and,  hence,  the  Raman  coherence  pbc  be¬ 
cause  of  population  transfers  to  |6')  via  RET  during  the 
long-delayed  period,  reducing  the  total  ERE-CARS  po¬ 
larization.  Hence,  the  threshold  intensity  for  saturation 
of  ERE-CARS  increases.  For  pump  and  Stokes  pulse  in¬ 
tensities  well  below  Raman  saturation  condition,  the  03- 
based  criterion  is  accurate  if  we  can  temporally  isolate 
the  determining  transition  of  interest  but  within  a  limit 
that  probe  delay  is  much  smaller  than  the  effective  decay 
rate  of  the  system.  However,  in  all  of  the  above  cases, 
the  ©3-based  criterion  predicts  the  saturation  threshold 
for  probe  intensity  accurately  within  one  order  of  magni¬ 
tude. 


V.  COMPARISON  OF  ERE-CARS 
SATURATION  IN  LONG-  AND  SHORT-PULSE 
DURATION  REGIMES 


In  the  previous  section,  we  demonstrated  in  Fig.  3  that 
the  threshold  intensity  of  saturation  of  Sere  is  a  func¬ 
tion  of  the  duration  of  the  probe  pulse  [as  also  shown  in 
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FIG.  5.  (Color  online)  (a)  Saturation  of  integrated  ERE- 
CARS  signal  intensity  for  different  pulse  durations  in  longer- 
pulse  regime  between  10  ns  and  10  ps.  (b)  Pulse  area  corre¬ 
sponding  to  peak  intensities  of  probe  pulse. 
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FIG.  6.  (Color  online)  (a)  Saturation  of  integrated  ERE- 
CARS  signal  intensity  for  different  pulse  durations  in  shorter- 
pulse  regime  between  10  ps  and  10  fs.  (b)  Pulse  area  corre¬ 
sponding  to  peak  intensities  of  probe  pulse. 


Eq.  (17)].  However,  from  the  earlier  literature  on  ERE- 
CARS,  we  know  that  if  the  peak  intensity  of  the  probe 
transition  is  increased  such  that  the  Rabi  frequency  of 
the  probe  is  larger  than  the  effective  decay  rate  of  the 
transition,  an  onset  of  a  Rabi  oscillation  between  the 
states  |  a)  o  |6)  engages  the  population  of  those  states 
through  the  interaction  with  the  probe  pulse.  If  we  in¬ 
crease  the  number  of  photons  further  (or  increase  the 
peak  intensity),  the  new  photons  do  not  have  an  oppor¬ 
tunity  to  interact  with  the  molecule  and,  hence,  can  no 
longer  be  absorbed,  leading  to  the  absorption  saturation 
in  | a)  1 6)  transition;  in  general,  this  is  also  referred 
as  the  transition  being  saturated.  For  CW  or  long-pulse 
excitations,  this  condition  can  be  achieved  if  the  driven 
system  reaches  an  underdamping  regime;  i.e. ,  the  peak 
Rabi  frequency  of  the  probe  pulse  is  larger  than  the  ef¬ 
fective  decay  rate  of  the  system,  i.e., 

^30  >  \/r ac<c  ,  (19) 

where,  rac  =  Ta  +  Tcc/  is  the  total  population  decay  in 
the  probe  transition.  It  is  clear  from  the  above  satura¬ 
tion  condition  that  in  the  long-pulse-excitation  regime 
the  threshold  intensity  of  the  saturation  is  independent 
of  the  duration  of  the  probe  pulse  T3,  unlike  in  the  ul¬ 
trafast  ERE-CARS  results  in  Fig.  3.  For  absorption  and 
fluorescence  with  short-pulse  excitation,  we  have  shown 
that  the  underdamping  condition  is  only  a  necessary  con¬ 
dition  [22].  The  ultrafast  saturation  of  absorption  and 


fluorescence  occurs  at  a  much  higher  intensity  when  the 
system  reaches  a  transient-gain  regime  within  the  dura¬ 
tion  of  the  driving  pulse.  Therefore,  for  ERE-CARS  it 
is  important  to  compare  the  saturation  of  the  signal  for 
different  probe-duration  regimes  to  investigate  what  are 
the  ranges  of  pulse  duration  where  the  saturation  cri¬ 
teria  transitions  from  a  purely  underdamping  condition 
for  long  pulses  to  the  transient-gain  condition  (03-based 
criteria)  for  the  ultrafast  saturation  of  the  ERE-CARS 
signal. 

In  the  following,  we  compare  the  saturation  of  the 
ERE-CARS  signal  over  an  extensively  wide  range  of  du¬ 
rations  of  the  probe  pulse  from  10  ns  to  10  fs.  For  the 
purpose  of  clarity,  we  have  plotted  Sere  for  73  —1  10  ns 
to  10  ps  in  Fig.  5(a)  and  t3  — >  10  ps  to  10  fs  in  Fig. 
6(a).  To  understand  the  validity  of  the  03-based  crite¬ 
ria  described  in  the  earlier  section,  we  have  also  plotted 
the  pulse  area  03  in  Figs.  5(b)  and  6(b),  correspond¬ 
ing  to  the  probe  intensities  used  in  Figs.  5(a)  and  6(a), 
respectively.  Note  that  in  Fig.  3,  we  maintained  the  du¬ 
rations  of  the  pump  and  Stokes  at  n  =  12  =  100  fs,  and 
only  the  probe-pulse  duration  r3  was  varied.  However, 
in  Figs.  5  and  6,  the  duration  of  all  three  pulses  is  con¬ 
sidered  to  be  the  same  T\  =  T2  =  t3  ,  and  they  are  varied 
simultaneously.  Furthermore,  real  ultrafast  ERE-CARS 
experiments  may  need  to  be  performed  with  a  probe  de¬ 
lay  to  avoid  any  non-resonant  background  in  the  observed 
signal;  but  to  ensure  uniformity  in  the  interaction  condi¬ 
tions  over  the  range  of  six-orders-of-magnitude  variation 


of  the  probe-pulse  duration,  we  maintain  a  zero  probe- 
pulse  delay  in  both  Figs.  5  and  6.  Also,  to  eliminate 
any  saturation  or  saturation-like  effect  in  Raman  coher¬ 
ence  caused  by  pump  or  Stokes  pulses  (as  in  Ref.  [21]), 
we  maintain  the  peak  intensities  of  the  pump  and  Stokes 
pulses  at  Iio  =  ho  =  1010  W/m2,  where  we  have  veri¬ 
fied  that  Raman  saturation  does  not  occur  even  in  the  ns 
regime. 

The  integrated  ERE-CARS  signal  Sere  for  long  pulses 
is  plotted  in  Fig.  5(a).  The  signal  increases  until 
ho  ~  1010  W/m2,  but  then  it  saturates  and  begins  to 
decrease  for  the  pulse  durations  from  10  ns  to  100  ps, 
which  corresponds  to  the  peak  Rabi  frequency  that  is 
larger  than  the  effective  decay  rate,  as  described  in  Eq. 
(19)  for  an  underdamping  regime.  Clearly,  the  thresh¬ 
old  intensity  for  saturation  is  independent  of  the  probe 
pulse  until  the  probe-pulse  duration  reaches  r3  ~  10  ps, 
unlike  the  dependence  described  in  Eq.  (17)  from  the  03 
criterion.  The  pulse  area  @3  corresponding  to  the  differ¬ 
ent  pulse  durations,  plotted  in  Fig.  5(b),  shows  that  03 
corresponding  to  the  pulse  durations  for  10  ns  and  100 
ps  are  ~  5007T  and  ~  507T,  respectively.  For  r3  =  100  ps 
and  t3  =  10  ps,  the  pulse  area  is  03  ~  4.57T  and  1.97T,  re¬ 
spectively.  Clearly,  for  long  pulses  with  durations  much 
longer  than  the  molecular  dephasing  timescales  (~  20 
ps),  the  pulse-area-based  criterion  completely  fails.  The 
saturation  threshold  is  determined  only  by  the  under¬ 
damping  condition  Eq.  (19).  As  noted  in  Refs.  [7,  16], 
a  significant  repumping  can  occur  during  the  duration 
of  the  excitation  pulse  that  has  been  really  helpful  for 
the  quenching  independence  of  the  ERE-CARS  signal  in 
NO.  However,  for  a  pulse  duration  on  the  order  of  the 
dephasing  timescale  of  the  molecule,  the  underdamping 
condition  is  not  a  sufficient  condition  for  saturation  of 
the  ERE-CARS  signal. 

We  have  plotted  the  ERE-CARS  signal  Sere  for  ultra- 
short  pulses  with  durations  ranging  from  10  ps  to  10  fs  in 
Fig.  6(a).  After  the  initial  rise  Sere  saturates  in  each 
of  the  cases  of  the  pulse  durations  of  the  probe  but  at 
very  different  ho,  unlike  those  in  Fig.  5.  The  signal  sat¬ 
uration  in  the  ultrafast  excitation  regime  with  the  probe 
pulse  duration  shorter  than  the  dephasing  time  scale  may 
be  understood  as  follows:  the  shorter  the  duration  of  the 
pulse  compared  to  the  effective  decay  timescale  of  the 
molecule,  the  higher  the  laser  intensity  needed  to  cause  a 
Rabi  cycle  to  occur  within  r3,  which  is  necessary  for  the 
saturation  of  the  probe  transition.  This  process  is  very 
similar  to  that  for  the  ultrafast  saturation  of  fluorescence 
described  in  Ref.  [22]  and  saturation  of  picosecond  po¬ 
larization  spectroscopy  [25,  26].  The  ultrafast  saturation 
of  S ere  follows  the  03-based  criterion  obtained  in  Eq. 
(17),  i.e. ,  the  threshold  intensity  of  saturation  of  Sere 
increases  quadratically  with  a  decrease  in  the  duration 
of  the  probe  pulse  r3.  The  pulse  area  corresponding  to 
each  r3  is  plotted  in  Fig.  6(b).  In  each  case  of  r3,  the 
saturation  of  S ere  follows  the  predicted  pattern  of  Eq. 
(17),  but  the  criterion  slightly  underpredicts  the  satura¬ 
tion  threshold  by  a  factor  of  two.  The  reason  is:  we  have 


assumed  the  probe  delay  to  be  zero  for  all  of  the  cases 
where  in  presence  of  all  three  fields,  some  active  repopu¬ 
lation  may  be  occurring  in  the  state  |  b)  and  |a),  pushing 
the  threshold  of  saturation  higher.  This  result  follows  the 
discussion  of  the  role  of  probe  delays  presented  in  the  pre¬ 
vious  section.  As  shown  earlier,  the  03-based  criterion 
will  be  exact  for  ultrafast  saturation  of  ERE-CARS  if 
the  probe  delay  is  maintained  within  two-to-three  times 
the  probe  pulse  duration.  The  predicted  threshold  of 
saturation  agrees  well  with  the  all-resonant  ERE-CARS 
experiments  carried  out  recently  [38]. 

It  may  be  noted  that  since  the  resonant  probe  pulse 
intensities  needed  for  saturation  are  very  high,  it  may 
cause  ac-Stark  shift  or  ionization  of  NO  via  two-photon 
transitions  from  X2n(u  =  1).  There  are  a  few  literature 
that  exist  on  the  ultrafast  UV-laser  excitation  of  NO,  e.g., 
Lopez-Marten  et  al.  have  shown  that  the  laser  intensi¬ 
ties  in  excess  of  1017  W/m2  with  140-fs  pulse  operating 
at  410  nm  can  cause  ac-Stark  shift  via  two-photon  tran¬ 
sition  [39],  (b)  laser  peak  intensity  of  2  x  1017  W/m2  at 
265  nm  (40-60  fs)  has  been  used  in  conjunction  with  a 
probe  at  400  nm  to  ionize  NO  via  resonant  excitation  of 
an  intermediate  level  to  observe  ac-Stark  shift  (see  Ref. 
[40]).  Though,  to  the  best  of  our  knowledge,  no  study  has 
been  reported  for  ionization  of  NO  by  UV  pulse  of  236 
nm  at  10  fs,  from  the  reported  experiments  of  NO  with 
200-fs  pump  at  800  nm  [41],  we  believe  no  significant 
ionization  would  not  occur  before  saturation  threshold 
of  ERE-CARS  signal  limit  is  reached,  even  for  the  10-fs 
excitation. 


VI.  SUMMARY 

In  summary,  we  obtained  a  simplified  criterion  for  the 
saturation  threshold  for  the  intensity  of  the  ultrafast  res¬ 
onant  probe  in  an  ERE-CARS  configuration.  Our  cri¬ 
terion  holds  rather  well  for  a  range  of  ultrashort  pulse 
durations  that  are  effectively  delayed  with  respect  to  the 
pump  and  Stokes  pulses.  We  have  compared  the  thresh¬ 
old  of  saturation  of  the  probe  between  pulse  durations 
ranging  from  10  ns  to  10  fs  and  shown  how  the  phys¬ 
ical  mechanism  of  saturation  differs  between  the  long- 
and  short-pulse  regimes  and,  hence,  differs  the  satura¬ 
tion  criterion.  The  underdamping  condition  is  a  suffi¬ 
cient  condition  for  saturation  of  ERE-CARS  in  the  long- 
pulse  regime,  leading  the  threshold  intensity  to  be  the 
same  for  up  to  100-ps  durations.  However,  for  ultrafast 
ERE-CARS,  an  additional  condition  must  be  satisfied 
that  the  probe  transition  must  changeover  from  tran¬ 
sient  absorption  to  transient  gain  regime;  hence,  with 
a  pulse-area-based  calculation,  we  obtained  a  criterion 
that  predictively  determines  the  ERE-CARS  saturation 
threshold  intensity  of  the  probe  pulse  for  durations  rang¬ 
ing  from  10  ps  to  10  fs  within  one  order  of  magnitude.  In 
the  ultrashort  pulse  excitation  regime,  the  saturation- 
threshold  intensity  is  quadratically  dependent  on  the 
inverse  of  the  pulse  duration,  unlike  in  the  long-pulse 
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regime.  The  pulse-area-based  study  may  be  generalized 
to  obtain  thresholds  for  saturation  of  multi-photon  pro¬ 
cesses  [42],  laser- induced  breakdown  [43],  ionization  or 
other  nonlinear  optical  processes  in  a  variety  of  systems 
that  employs  the  ultrafast  excitation,  and  also  for  differ¬ 
ent  pulse  shapes  [44,  45]. 


D-2329  and  FA8650-15-D-2518)  and  the  U.S.  Air  Force 
Office  of  Scientific  Research  (Dr.  Enrique  Parra,  Pro¬ 
gram  Manager).  This  manuscript  has  been  cleared  for 
public  release  (No.  88ABW-2015-3756). 


VII.  ACKNOWLEDGMENT 

Funding  for  this  study  was  provided  by  the  U.S.  Air 
Force  Research  Laboratory  (Contract  Nos.  F33615-03- 


[1]  S.  Roy,  J.  R.  Gord,  and  A.  K.  Patnaik,  Prog.  Energy 
Combust.  Sci.  36,  280  (2010). 

[2]  S.  F.  Hanna,  W.  D.  Kulatilaka,  Z.  Arp,  T.  Opatrny,  M. 
O.  Scully,  .1.  P.  Kuehner,  and  R.  P.  Lucht,  Appl.  Phys. 
Lett.  83,  1887  (2003). 

[3]  S.  A.  J.  Druet,  B.  Attal,  T.  K.  Gustafson,  and  J.  P. 
Taran,  Phys.  Rev.  A  18,  1529  (1978). 

[4]  B.  Attal-Tretout,  S.  C.  Schmidt,  E.  Crete,  P.  Dumas,  and 

J.  P.  Taran,  J.  Quant.  Spectrosc.  Radiat.  Transf.  43,  351 
(1990). 

[5]  M.  Schenk,  T.  Seeger,  and  A.  Leipertz,  Appl.  Opt.  44, 
4157  (2005). 

[6]  S.  Roy,  W.  D.  Kulatilaka,  S.  V.  Naik,  N.  M.  Laurendeau, 
R.  P.  Lucht,  and  J.  R.  Gord,  Appl.  Phys.  Lett.  89,  104105 
(2006). 

[7]  A.  K.  Patnaik,  S.  Roy,  J.  R.  Gord,  R.  P.  Lucht,  and  T. 
B.  Settersten,  J.  Chenr.  Phys.  130,  214304  (2009). 

[8]  W.  D.  Kulatilaka,  N.  Chai,  S.  V.  Naik,  N.M.  Laurendeau, 
R.  P.  Lucht,  J.  P.  Kuehner,  S.  Roy,  and  J.  R.  Gord,  Opt. 
Lett.  31,  3357  (2006). 

[9]  N.  Chai,  W.  D.  Kulatilaka,  S.  V.  Naik,  N.  M.  Lauren¬ 
deau,  R.  P.  Lucht,  J.  P.  Kuehner,  S.  Roy,  V.  R.  Katta, 
and  J.  R.  Gord,  Appl.  Phys.  B:  Lasers  Opt.  88,  141 
(2007). 

[10]  W.  D.  Kulatilaka,  N.  Chai,  S.  V.  Naik,  S.  Roy,  N.  M. 
Laurendeau,  R.  P.  Lucht,  J.  P.  Kuehner,  and  J.  R.  Gord, 
Opt.  Cornmun.  274,  441  (2007). 

[11]  N.  Chai,  S.  V.  Naik,  N.  M.  Laurendeau,  R.  P.  Lucht,  S. 
Roy,  and  J.  R.  Gord,  Appl.  Phys.  Lett.  93,  091115  (2008) 

[12]  J.  P.  Kuehner,  S.  V.  Naik,  W.  D.  Kulatilaka,  N.  Chai, 
N.  M.  Laurendeau,  R.  P.  Lucht,  M.  O.  Scully,  S.  Roy,  A. 

K.  Patnaik,  and  J.  R.  Gord,  J.  Chem.  Phys.  128,  174308 
(2008). 

[13]  A.  C.  Eckbreth,  Laser  Diagnostics  for  Combustion  Tem¬ 
perature  and  Species  (Gordon  and  Breach,  Amsterdam, 
1996). 

[14]  N.  Chai,  R.  P.  Lucht,  W.  D.  Kulatilaka,  S.  Roy,  and  J. 
R.  Gord,  J.  Chem.  Phys.  133,  084310  (2010). 

[15]  L.  A.  Rahn,  R.  L.  Farrow,  M.  L.  Koszykowski,  and  P.  L. 
Mattern,  Phys.  Rev.  Lett.  45,  620  (1980). 

[16]  A.  K.  Patnaik,  S.  Roy,  R.  P.  Lucht,  and  J.  R.  Gord,  J. 
Mod.  Opt.  55,  3263  (2008). 

[171  R.  Altkorn  and  R.  Zare,  Annu.  Rev.  Phys.  Chem.  35, 
265  (1984). 

[18]  R.  P.  Lucht,  R.  L.  Farrow,  and  D.  J.  Rakestraw,  J.  Opt. 
Soc.  Am.  B  10,  1508  (1993). 


[19]  O.  Kinrot  and  Y.  Prior,  Phys.  Rev.  A  50  (3),  R1999 
(1994) 

[20]  R.  P.  Lucht  and  R.  L.  Farrow,  J.  Opt.  Soc.  Am.  B  5, 
1243  (1988). 

[21]  A.  K.  Patnaik,  S.  Roy,  and  J.  R.  Gord,  Phys.  Rev.  A  87, 
043801  (2013). 

[22]  A.  K.  Patnaik,  S.  Roy,  and  J.  R.  Gord,  Phys.  Rev.  A  90, 
063813  (2014). 

[23]  R.  W.  Boyd,  Nonlinear  Optics.  (Academic  Press,  San 
Diego,  CA,  1992). 

[24]  W.  Demtroder,  Laser  Spectroscopy,  Vol.  2:  Experimental 
Techniques.  (Springer- Verlag,  Berlin/Heidelberg,  2008). 

[25]  T.  A.  Reichardt,  F.  D.  Teodoro,  R.  L.  Farrow,  S.  Roy, 
and  R.  P.  Lucht,  J.  Chem.  Phys.  113,  2263  (2000). 

[26]  Sukesh  Roy,  Robert  P.  Lucht,  and  Thomas  A.  Reichardt, 
J.  Chem.  Phys.  116,  571  (2002). 

[27]  S.  L.  McCall  and  E.  L.  Hahn,  Phys.  Rev.  Lett.  18,  908 
(1967). 

[28]  G.  M.  Gen  kin,  Phys.  Rev.  A  58,  758  (1998). 

[29]  R.  Kienberger  and  F.  Krausz,  ’’Subfemtosecond  XUV 
pulses:  Attosecond  metrology  and  spectroscopy”  in  Few- 
Cycle  Laser  Pulse  Generation  and  Its  Applications, 
Ed.  F.  X.  Kartner  (Springer- Verlag,  Bcrlin/Heidelberg, 
2010). 

[30]  L.  Yuan,  G.  O.  Ariunbold,  R.  K.  Murawski,  D.  Pestov, 

X. Wang,  A.  K.  Patnaik,  V.  A.  Sautenkov,  A.  V.  Sokolov, 

Y.  V.  Rostovtsev,  and  M.  O.  Scully,  Phys.  Rev.  A  81, 
053405  (2010). 

[31]  G.  Heck,  A.  Filin,  D.  A.  Romanov,  and  R.  J.  Levis,  Phys. 
Rev.  A  87,  023419  (2013). 

[32]  Y.  Yonemaru,  A.  F.  Palonpon,  S.  Kawano,  N.  I.  Smith, 
S.  Kawata,  and  K.  Fujita,  Phys.  Rev.  Appl  4,  014010 
(2015). 

[33]  G.  B.  Arfken  and  H.  J.  Weber,  Mathematical  Methods  for 
Physicists  (Sixth  Edition,  Academic  Press,  Amsterdam, 
2005). 

[34]  M.  Abramowitz  and  I.  A.  Stegun,  Handbook  of  Meth- 
matical  Functions  (Dover  Publications,  Inc,  1972). 

[35]  B.  D.  Prince,  A.  Chakraborty,  B.  M.  Prince,  and  H.  U. 
Stauffer,  J.  Chem.  Phys.  125,  044502  (2006). 

[36]  D.  Pestov,  R.  K.  Murawski,  G.  O.  Ariunbold,  X.  Wang, 
M.  Zhi,  A.  V.  Sokolov,  V.  A.  Sautenkov,  Y.  V.  Rostovt¬ 
sev,  A.  Dogariu,  Y.  Huang,  and  M.  O.  Scully,  Science 
316,  265  (2007). 

[37]  H.  U.  Stauffer,  J.  D.  Miller,  M.  N.  Slipchenko,  T.  R. 
Meyer,  B.  D.  Prince,  S.  Roy,  and  J.  R.  Gord,  J.  Chem. 


10 


Phys.  140,  024316  (2014) 

[38]  S.  Roy,  H.  U.  Stauffer,  J.  Schmidt,  P.  Wrzesinski,  and  J. 
R.  Gord  (to  be  published). 

[39]  R.  B.  Lopez-Martens,  T.  W.  Schmidt,  and  G.  Roberts, 
Phys.  Rev.  A  62,  013414  (2000). 

[401  Z.  Sun,  H.  Liu,  N.  Lou,  and  S.  Cong,  Chem  Phys.  Lett. 
369,  374  (2003). 


[41]  A.  Talebpour,  S.  Larochelle,  and  S.  L.  Chin,  J.  Phys.  B: 
At.  Mol.  Opt.  Phys.  30,  1927  (1997). 

[42]  H.  U.  Stauffer,  W.  D.  Kulatilaka,  J.  R.  Gord,  and  S.  Roy, 
Opt.  Lett.  36,  1776  (2011). 

[43]  E.  W.  Van  Stryland,  M.  J.  Soileau,  A.  L.  Smirl,  and  W. 
E.  Williams,  Phys.  Rev.  B.  23,  2144  (1981). 

[44]  P.  R.  Berman,  L.  Yan,  K.  -H.  Chiam,  and  R.  Sung,  Phys. 
Rev.  A  57,  79  (1998). 

[45]  A.  K.  Patnaik,  S.  Roy,  and  J.  R.  Gord,  Phys.  Scr.  T165 
014031  (2015). 


